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Abstract 

A dedicated analysis of the muon-induced background in the EDELWEISS dark matter search has been 
performed on a data set acquired in 2009 and 2010. The total muon flux underground in the Laboratoire 
Souterrain de Modane (LSM) was measured to be — (5.4 ± 0.2tQ g) muons/m^/d. The modular design 
of the /i-veto system allows the reconstruction of the muon trajectory and hence the determination of the 
angular dependent muon flux in LSM. The results are in good agreement with both MC simulations and 
earlier measurements. Synchronization of the /i-veto system with the phonon and ionization signals of the 
Ge detector array allowed identification of muon-induced events. Rates for all muon-induced events F'' = 
(0.172 ±0.012) evts/(kg • d) and of WIMP-like events F^"" = 0.008tg !][J| evts/(kg • d) were extracted. After 
vetoing, the remaining rate of accepted muon-induced neutrons in the EDELWEISS-II dark matter search 
was determined to be F[J.~^ < 6 • 10~^evts/(kg • d) at 90%C.L. Based on these results, the muon-induced 
background expectation for an anticipated exposure of 3000 kg • d for EDELWEISS-III is A'^^ooofcg d < ^-^ 
events. 
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1 Introduction 

Astrophysical and cosmological observations i38l [29] 
have led to the current cosmological concordance 
model that incorporates both dark energy and dark 
matter as driving forces in the evolution of the Uni- 
verse. Dark matter herein constitutes a matter-like 
fluid ruling the dynamics of structure formation as 
well as significantly affecting the dynamics of our own 
galaxy today [27j. Assuming extensions of the Stan- 
dard Model with a new particle originally in ther- 
modynamic equilibrium, but decoupling shortly af- 
ter the Big Bang, a promising class of dark matter 
candidates called WIMPs (Weakly Interacting Mas- 
sive Particles) arises. An annihilation cross section of 
the order of the weak scale is needed to produce the 
observed cosmic abundance [18], rather independent 
of the WIMP mass which is typically assumed to be 
of 0(100 GeV). With these properties, this candidate 
is non-relativistic at freeze-out and hence consistent 
with our understanding of structure formation. 

The challenge of direct dark matter search experi- 
ments is the detection of the extremely rare WIMP- 
nucleus scattering events. Furthermore the expected 
energy spectrum from WIMP-nucleus scattering ex- 
hibits an exponential distribution, tailing off towards 
higher energy and can be easily mimicked by back- 
ground sources. Thus, in order to ensure that a signal 
is caused by WIMPs one has to understand all back- 
ground sources and to suppress as many as possible. 
The overwhelming majority of observed events are 
gamma and beta particles from ambient radioactiv- 
ity. These are rejected with the powerful discrimi- 
nation capabilities offered by the simultaneous read- 
out of two channels such as phonons and scintillation 
[3 O |20] , or phonons and ionization [H [3] , or scintilla- 
tion and ionization [391 13 E] > ^^'^ fiducialization (bulk 
- surface) in some detectors. In contrast, a neutron 
background causes a WIMP-like signal (a nuclear re- 
coil) in the detectors and can only be discriminated 
via multiple scattering of neutrons or coincident fi- 
veto signals if induced by cosmic muons. 

At levels of sensitivity to spin-independent elas- 
tic WIMP-nucleon cross sections of a few 10"^'* cm^ 



(CDMS and EDELWEISS with 614 kg • d of exposure 
[3]) or lower (XENONIOO [9] with 2323.7 kg • d equiv- 
alent exposure), further reducing the background is 
a demanding task. It requires a careful analysis of all 
known backgrounds and the extrapolation to the next 
stages of experiments [H]. For EDELWEISS, promi- 
nent sources of backgrounds are a potential failure 
to separate electron recoils from nuclear recoils via 
the charge and phonon measurement, and neutrons 
produced via {a, n) reactions, spallation as well as 
neutrons induced by muons. Old results have con- 
siderable systematic uncertainties due to the lack of 
detailed Monte Carlo modeling of muon transport, 
secondary particle production and detection within a 
specific detector configuration. More recent experi- 
ments (see for instance [H [lOl |T6l [35]) provide more 
reliable data for different targets but some of them 
are still not consistent with the existing models based 
on the GEANT4 toolkit [2] or FLUKA code ^ (see 
also simulations and discussions in Refs. [TT| \M\). 

In this article we present the results of an investiga- 
tion of the muon-induced background in the EDEL- 
WEISS experiment. We start with an introduction 
into the experimental setup, describing briefly the 
Ge bolometers and focussing on the fi-veto system 
and its data acquisition (Section [5]). In Section [31 
we derive the muon detection efficiency and recon- 
struct muon information such as angular distribu- 
tions and track coordinates. Based on these data, 
a search for coincidences of muon veto and bolome- 
ter events has been performed with results presented 
in Section [H We conclude with an outlook for the 
muon-induced background expectation for the up- 
coming EDELWEISS-III dark matter search. 

2 Experimental setup 

To shield against cosmic backgrounds, the EDEL- 
WEISS experiment is housed in the Laboratoire 
Souterrain de Modane (LSM), an underground lab- 
oratory under the French Alps with an average rock 
overburden of 4800 m w. e. (meters water equiv- 
alent). This attenuates the muon flux by six or- 
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ders of magnitude compared to sea level down to 
~ 5/x/m^/d. The experimental setup depicted in 
Fig. [1] consists of several layers of active and passive 
shielding surrounding the innermost part, a dilution 
cryostat. This cryostat with a working temperature 
of 18 mK contains the ultrapure Ge crystals used as 
detectors for the WIMP search. 

The outermost shielding is an active /i-veto sys- 
tem of plastic scintillator modules, described in Sec- 
tion [221 It is mounted on a 50 cm thick polyethylene 
(PE) layer, inside which the actual cryostat hous- 
ing is shielded with 20 cm thick lead, of which the 
inner 2 cm are ancient lead (Pb). The PE moder- 
ates the neutron energies to sufficiently low levels for 
a significant fraction of them to be absorbed. The 
innermost 2 cm of the lead were salvaged from the 
wreck of an ancient Roman galley |32) in order to 
reduce the 46.5 keV 7-line of the ^^'^Pb as well as 
X-rays and bremsstrahlung photons from the /3~ de- 
cay of ■^^'^Bi. The Ge bolometers are described in 
Section 12.11 Their signals are sampled using a com- 
mon clock with a 100 kHz frequency. That same clock 
is transferred to the ^-veto system data acquisition 
(DAQ) via optical fibers. Apart from this synchro- 
nization, both systems run individually and muon- 
induced bolometer events are identified offline within 
the main analysis p!3]. An additional 1 ton liquid 
scintillator detector designed to study muon-induced 
neutrons with better statistics was installed at the 
end of 2008 (Fig. [1^, detector outside the shielding 
at the lower right). The detailed setup and first re- 
sults obtained with this additional detector are de- 
scribed in ['30j. Further measurements of the ambient 
neutron flux are performed inside and outside of the 
shielding with bare ^He neutron counters j42[ I23j. 

2.1 Bolometer array and DAQ 

The Ge crystals for EDELWEISS-II are grown from 
ultrapure Ge with less than 10^^ impurities per cm^ 
and cut as cylinders of 70 mm diameter and 20 mm 
height. The simultaneous readout of a heat signal 
with a neutron transmutation doped (NTD) Ge sen- 
sor [37j l36j ) and the ionization with charge collection 



electrodes enables the distinction between electron 
and nuclear recoils. During the data taking period 
2009-2010, the dilution cryostat provided a stable de- 
tector temperature of 18 mK for more than one year. 
Five ID detectors with beveled edges and a mass of 
360 g each were used together with flve plain cylin- 
drical detectors of 410 g mass. These ID detectors 
were equipped with interdigitized collecting and veto 
electrodes [TU]- The Al electrodes are evaporated on 
the top and bottom of the Ge crystal. Each ring 
has a width of 200 fiui and 1.8 mm pitch, and is typ- 
ically biased at voltages of ±4V for the collecting 
and =f1.5V for the veto electrode. Solid guard elec- 
trodes covering the outer rim and the beveled faces 
allow suppression of events in the outer regions of 
the crystals. The electrical fleld lines run orthogonal 
to the surface for bulk events, but parallel for surface 
events. This design therefore allows to exclude events 
in the surface layers of the detector, which are prone 
to charge collection deficiencies. 

In the analysis the discrimination between electron 
and nuclear recoils is based on the ionization yield Q, 
which is the ionization energy (calibrated for electron 
recoils with the 356 keV line of ^^^Ba) relative to the 
recoil energy. For electron recoils the value of Q is 
set to unity by calibration. Nuclear recoils trans- 
mit most of their energy directly to vibrations of the 
crystal lattice (phonons), without going through ion- 
ization processes. Thus, they generally have a much 
smaller ionization yield than electron recoils. For the 
EDELWEISS detectors, the ionization yield of nu- 
clear recoils as a function of energy follows closely 
the prediction of the Lindhard theory [33] and can be 
parameterized through the expression, determined in 
an extensive neutron calibration |37[ \^3\ : 

QLin = 0.16-(ERec[keV])0-i8. (I) 

Assuming a Gaussian uncertainty of the Q value dis- 
tribution with energy-, detector- and time-dependent 
resolutions fTion(Eion, t, det) and (Theat(Eheat, t, det), a 
recoil band can be defined based on the standard de- 
viation 

C'Q = 0-Q(£^Rec,(Tion(t,det),CTheat(t,det)). (2) 
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(a) 



(b) 



Figure 1: (a) Cross-section of the EDELWEISS experimental setup as implemented in the MC simulations, (b) Schematic view 
of the modular /i-veto system. 



For the ID detectors, rejection powers of better than 
1 in 10^ for 7 rays in the bulk and 1 in 6-10^ for 
surface events were determined using a ^^^Ba 7 and 
a ^"^'^Pb /3 source, respectively [19]. Whilst these re- 
jection powers are very good, the surface event re- 
jection together with the rejection of events on the 
radial surface, so-called guard events, significantly 
reduces the volume of accepted events. The bulk 
or fiducial volume, where reliable charge collection 
is expected, was defined in the following way. Each 
of the signals on veto and guard electrodes, as well 
as the charge imbalance between the two collecting 
electrodes had to be consistent with zero at 99 %C.L., 
according to the resolution of these parameters. The 
efficiency of those cuts was measured using the low- 
energy 8.98 keV and 10.34 keV peaks associated with 
the cosmogenic activation of germanium, producing 
a uniform contamination with the isotopes ^^Zn and 
^^Ge. Averaged over time and all detectors, a typical 
fiducial mass of (166 ± 6) g per detector was found 



With a live data taking period of 307 days dur- 
ing the cryogenic run from March 2009 to May 2010, 
an effective fiducial exposure of (481 it 17) kg • d has 
been collected for the analysis of the muon-induced 
background in the bolometers including a factor of 
0.9 to account for the region of interest (Rol) being 
the 90%C.L. nuclear recoil band. This exposure is 
larger than the 384 kg • d used in the WIMP search 
in [13], as the present study of coincidences of muon 
veto and germanium detector events does not require 
the same strict quality cuts on the noise levels in the 
germanium detectors. While the fiducial volume cut 
is paramount for the DM search |13j it can be bene- 
ficial to include less well reconstructed events closer 
to the surface of the detectors to increase statistics 
for coincidence studies. Therefore, we also present 
the results of a coincidence study using looser cuts 
for events in the whole volume of the detectors, for 
which the exposure amounts to 1504 kg-d. 
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2.2 Muon veto system 

The outermost layer of the EDELWEISS setup is the 
fi-veto system, a layer of plastic scintillator acting 
as an active shield to tag any muon-induced activ- 
ity in the bolometer system. It consists of 46 indi- 
vidual plastic scintillator modules (Bicron BC-412) 
of 65 cm width, 5 cm thickness and lengths of 2 m, 
3.15 m, 3.75 m and 4 m adding up to a total surface 
of 100 m^. The modules are assembled in a stainless 
steel mounting frame and attached to the PE shield- 
ing to ensure an almost hermetic coverage of the ex- 
periment setup. Due to the need for access to the 
cryostat during maintenance intervals, parts of the fi- 
veto system are positioned on rails on top of the two 
movable PE blocks to allow a roll-back during cryo- 
stat access (Fig. [1]). Whereas the upper level (Nl) 
of the //-veto system almost completely covers the 
PE shielding, except for small openings for rails, the 
part below the cryostat level (NO) has more promi- 
nent gaps due to the cryogenic supply lines, and the 
pillars on which the entire experiment is mounted. 
The overall geometrical muon flux coverage of the 
//-veto system amounts to 98%. To ensure proper 
closure during data taking, the exact positions of the 
two mobile veto parts are monitored via regular po- 
sition measurements using laser interferometry. 

Each module is viewed by two groups of four 2-inch 
PMTs (Philips Valvo XP2262/PA). Each group has 
an individual high voltage (HV) supply and readout, 
leading to a total of 92 HV and signal channels. Sim- 
ilar modules have been used as a muon veto counter 
in the KARMEN experiment [12]. A complete de- 
scription of the modules including measurements of 
the effective attenuation length and spectral quan- 
tum efficiency of the scintillator can be found in [40j . 
The data acquisition is based on VME electronics 
with special cards developed in- house [25]. An event 
recording in the veto system is triggered once the 
two signals of a module each pass a trigger threshold 
within a coincidence window of 100 ns. With leading 
edge discriminators set at the same level of 150 mV, 
individual modular thresholds are achieved by adjust- 
ing the HV values. Within an event, all hits above 



threshold in the modules are then registered within 
an interval of ±125 ns around the trigger time. For 
each signal channel, the individual time of a hit is 
recorded in a 128-channel TDC (CAEN V767) and 
the analog signals are converted in three 32-channel 
QDCs (CAEN V792). After the coincidence window, 
any further input is blocked for about 50 fis [25j until 
the event data are read out via the VME bus and all 
active elements are reset by the acquisition software. 
The time information with respect to the bolome- 
ters is received via an optical fiber which transmits 
the 10 IIS time stamp in 48 bits of 64 ns width each. 
This sequence of bits is converted via an FPGA-based 
board and then read out via the VME bus. To moni- 
tor the overall status of the veto modules, the hit rate 
of all PMT groups (not requiring any coincidence) is 
recorded in scalers. These data are processed contin- 
uously to check the performance of the /i-veto system 
via a graphical web interface. 

A muon traversing a scintillator module as a mini- 
mum ionizing particle deposits on average more than 
lOMeV (see Fig. [2] for a typical Landau spectrum). 
The average energy deposition is 11.8 MeV for hor- 
izontal modules on the top and bottom, and up to 
24 MeV for modules on the sides [26]. As seen in 
Fig. [21 these values are above typical energy deposi- 
tion from the ambient radioactivity. For each chan- 
nel, the energy deposition is stored as the integrated 
charge Q in units of ADC values. With the compar- 
ison of the Landau most probable values (MPV) of 
the integrated charge Q of specific two module coinci- 
dence cuts with dedicated MC simulations, a module- 
individual calibration is obtained with a mean cali- 
bration coefficient of 5.44 keV/ ADC- value. The ef- 
fective threshold of the muon veto modules has been 
chosen to be much lower than the MPVs (about 
5 MeV at the center of a module, rising towards the 
module ends due to the position-dependent light out- 
put) to be also able to detect grazing muons and thus 
to keep the muon detection efficiency as high as pos- 
sible. As a result, the rate of recorded events in the 
entire system is of the order of one per second, to be 
compared to an expected rate of throughgoing muons 
hitting at least one module per system section Nl and 



5 




10 20 30 

Energy (MeV) 



Figure 2: Distribution of energy depositions in a top module. 
The histogram (solid black) shows events recorded for a stan- 
dard trigger condition (internal coincidence of the module), 
scaled by 0.01. Data points (in red) describe energy deposi- 
tions of muon candidates selected by requiring full non-zero 
ADC and TDC information in another module (two modules 
in coincidence) and follow a Landau distribution. 

NO of ~ 3.5 X 10~^ Hz. Whereas the rate of recorded 
events is dominated by ambient background, candi- 
dates of throughgoing muons can be easily extracted 
from the data by requiring two modules in coinci- 
dence. 

Vetoing a bolometer event due to activity in the 
fi-veto system is performed in the offline analysis. 
A veto is applied in a conservative time window of 
±1 ms around each event recorded in the //-veto sys- 
tem, regardless of the fact that these events are dom- 
inated by ambient background. This offline reduc- 
tion of effective exposure due to the entire /i-veto 
system results in a deadtime of the order of 0.2 %. 
As the veto system is completely independent of the 
bolometer system and its data acquisition, data of 
these scintillator modules are acquired continuously, 
only interrupted by maintenance intervals. 

For each event in the /i-veto system, the position of 
an energy deposition can be reconstructed using the 
module number and the time difference between the 
signal arrival at each module end. The time is stored 
in TDC units, where one TDC unit corresponds to 
0.8 ns. A typical distribution of the time difference 
within a ^-veto module is shown in Fig. [3l The 
black distribution recorded with the standard trig- 
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Figure 3: Time difference of PMT signals within a top mod- 
ule for events with standard trigger condition (internal coinci- 
dence) in solid black, scaled by 0.01, and for muon candidates 
(red data points), selected with a two module coincidence cut. 

ger condition shows the threshold dependence along 
the module length with the minimum threshold and 
largest count rate at the center of the module. Ad- 
ditional contributions are visible at At ~ ±30 ns due 
to a thicker scintillator and thus higher detection ef- 
ficiency at the module ends. The tails to much larger 
time differences are due to the limited time resolution 
for low energy deposition as well as to light emission 
in the PMTs. In contrast, muon candidates show a 
much more confined distribution of time differences 
as can be seen in Fig. [3l The rather flat distribu- 
tion with sharp edges reflects the physical length of a 
module with homogenously distributed hits along the 
module axis. Thus, with the event data on energy de- 
position and time, sufficient information is provided 
to perform cuts selecting muon candidates and to re- 
construct muon trajectories (see Section [33]) . 

3 Detecting muons 

In this section, we derive the detection efficiency of 
the //-veto system in two different ways. Based on a 
detailed simulation of muons and their interactions 
in the EDELWEISS setup (Section 13. 1|) using indi- 
vidual modular efficiency curves determined in situ 
from data (Section 13. 2p . the first method is a combi- 
nation of MC results and data (Section 13. 3p . The sec- 
ond method uses bolometer data identified as muon- 
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Figure 4: LSM geometry and muon generation as implemented 
in Geant4. The hemisphere for muon generation is placed in 
the rock surrounding the laboratory hall with the EDELWEISS 
and NEMO setups. The NEMO-III experiment is implemented 
as simplified PE shields whereas the EDELWEISS-II setup is 
implemented with all details (see also Fig. [TJ. 

induced events searching for a coincident signal in 
the //-veto system (Section I3.4p and thus is entirely 
based on data extracted during the WIMP search 
itself. The measured angular distribution of muons 
(Section 13. 5p as well as the total muon flux are de- 
termined (Section l3.6p . 

3.1 Simulation of muons in LSM 

Modeling of muons and their interactions is per- 
formed with the Geant4 package version 4.9.2 OIU]. 
For these simulations, the detailed geometry of the 
EDELWEISS setup as given in Fig. [T] has been ex- 
tended to incorporate nearby structures in the LSM 
laboratory, namely the rock surrounding the labora- 
tory and the adjacent NEMO-III experiment. Fig- 
ure H] shows the implemented geometry as well as the 
hemisphere at a distance of 30 m from the center of 
the experiment on which muons are generated. The 
full MC simulation consists of 2- 10^ generated muons 
in the energy range from 2 GeV-200 TeV. 40.9% of 
the started muons have a momentum vector that ge- 
ometrically would hit the /i-veto system, the rest illu- 
minates the close surroundings. However, only 33.9% 
of the muon tracks actually enter the veto volume. 
The difference of 7% is due to muons stopped within 



the rock overburden. 

The energy and angular dependent muon flux 
above the laboratory is generated according to the 
energy dependent muon spectrum at sea level p^ : 

dNp _ 0.14(^^/GeV)'^ 
dEo cm^ s sr GeV 



I , I i.iEf,cose , i.ig^cose I ' ^ 

V + 115 GeV ^ + 850 GeV / 

with the cosmic ray exponent 7 = 2.7, which is then 
modified due to energy loss of muons and the shield- 
ing in the rock overburden of LSM. This rock over- 
burden h{@, <I>) is taken from a detailed elevation map 
with 1° X 1° resolution 

The muons are generated according to the result- 
ing angular dependent muon flux on a half sphere 
of 30 m radius centered around the cryostat of the 
EDELWEISS-II experiment. In order to ensure a 
homogeneous illumination of the /i-veto system, the 
starting point of a muon with momentum vector p 
is randomized within a distance 5 < 5 m of a ref- 
erence point on the hemisphere. After generation, 
particles are propagated through the entire setup 
based on the QGSP_BIC_HP reference physics list 
of Geant4 version 4.9.2 with an adaptation suggested 
by [26]. This includes the use of the "low energy ex- 
tension" packages instead of the standard ones for 
electron, positron and gamma interactions. Further, 
G4MuNuclear was added to model muon-hadronic 
interactions, and a CHiral Invariant Phase Space 
(CHIPS) based model was used for gamma-nuclear 
physics. 

To model hadronic interactions, especially inelas- 
tic neutron scattering, two further changes were ap- 
plied. At energies between 20MeV and 70MeV, the 
package G-^PreCompoundModel (PC) mediates the 
transition from G4NeutronHPInelastic (HP, up to 
20MeV) to G4BmaryCascade (BIG, from 70MeV to 
6 GeV). To bridge a model gap at energies 6 GeV < 
En.p ^ 12 GeV between BIG and the high energy 
Quark Gluon String Precompound (QGSP) model 
implemented in G4QGSModel, the Low Energy Pa- 
rameterized (LEP) model G4LENeutronInelastic was 
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inserted. This description of the physical processes 
has been proven to be effective in terms of CPU con- 
sumption and to be accurate when comparing with 
experimental results |26l [28] . Using these interaction 
models, the distance of 30 m of rock between starting 
point and experimental setup allows the development 
of hadronic and electromagnetic showers into equilib- 
rium and ensures the proper abundance of secondary 
particles in the simulation. Energy depositions by 
the primary muon and by all secondary particles are 
then stored for each individual veto module as well as 
for each bolometer in the cryostat. To avoid bound- 
ary effects, two additional meters of rock below the 
entire setup are included in the geometry. 

3.2 Trigger efficiency of individual veto 
modules 

To determine the detection efficiency for a given sim- 
ulated energy deposition in a module, the individual 
trigger threshold of all modules must be known. Al- 
though this is a position dependent function along 
the axis for a given module i, it can be measured in 
situ as an efficiency ei{E) averaging over all hit posi- 
tions. Such a measurement provides an accurate effi- 
ciency curve accounting for the local environment of 
the underground laboratory, ageing of detector com- 
ponents over the long term measurement as well as 
adjustments of the high voltage applied to the PMTs. 

The measurement of ei{E) makes use of the fact 
that, having an internal coincidence in an individual 
module (TDC data from both ends within a 100 ns 
window), all non-zero ADC and TDC data of other 
modules are also read out and stored. This additional 
information allows to deduce the trigger efficiency in 
the following way: We select all events where there is 
an energy deposition in a certain module recorded as 
non-zero ADC value(s), but where the trigger condi- 
tion was fulfilled by another module. Whenever both 
TDC channels of the module under consideration are 
non-zero, this module would have triggered as well. 
Otherwise, despite some energy deposition recorded 
in the ADC, the signals would not have triggered the 
event recording. The fraction of potential triggers 
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Figure 5: Example of the trigger efficiency ei (E) as a function 
of the energy deposition E in MeV in a module i (shown for a 
top module, with binomial error estimates). The green band 
is the la region of the efficiency fit, fixing the efficiency level 
for large energies to ^i{E) — 1. The functional form is given 
by the convolution of a threshold step- function with a gaussian 
energy resolution. 

over all events function of the calibrated energy 
is shown in Fig. [5] as an example of the trigger effi- 
ciency €i{E) of a specific module i. 

For muons, one can deduce a detection efficiency ef 
by weighting the energy dependent trigger efficiency 
with the measured Landau spectrum of each module: 

IZevN,dE)dE ■ ^> 

At energies below 6 MeV, the selection criterion for 
a muon candidate, i.e. two non-adjacent modules in 
coincidence, may also select low-energy secondaries in 
one of the modules instead of throughgoing muons. 
The energy cut of 6 MeV was chosen since the Landau 
spectra have a negligible fraction of events below that 
value in contrast to the data with a significant part 
due to secondaries. Fig. [6] shows the Landau spectra 
before and after correction of the efficiency, N^^i(E) 
and N;i^-{E) = N^4E)/ei{E), respectively. 

The accuracy of the individual muon detection effi- 
ciency is dominated by the uncertainty of the energy 
calibration based on fitting Landau spectra. The 
expected Landau MPV varies for different modules 
due to the mountain profile, the orientation of the 
modules and the muon selection cut and amounts to 
(-E-MPv) = 16 MeV. The Landau distribution itself 
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Figure 6: Measured energy deposition A^,t,i (E) of muon candi- 
dates for a top module requiring a second module with TDC 
coincidence (black); Energy spectrum N^^i{E) after correcting 
for inefficiencies (red). Both distributions have been fit with a 
Landau distribution excluding events at low energy which are 
due to energy depositions of secondary particles. 

is significantly broadened in comparison to that of 
vertically throughgoing muons. The energy calibra- 
tion per module is based on a modular comparison of 
MPV's from data fits (in ADC units) with MC sim- 
ulations in MeV and is estimated to be accurate at 
the 20% level. The resulting distribution of energies 
where the individual module trigger efficiency reaches 
50% is Gaussian with {E{ei = 0.5)) = 7 MeV as mean 
and a standard deviation of (t{E) = 1.9 MeV. Aver- 
aging the modular detection efficiency in Eq. U] 
results in (ef) = 95% with o-(ef ) = 4%. 

3.3 Muon detection efficiency derived 
from MC simulations 

In order to evaluate the muon detection efficiency 
for the entire //-veto system we defined two reference 
volumes for entering muons. At first, we take the 
volume defined by the //-veto system itself. A simu- 
lated energy deposition ii^ in a module i leads, with a 
probability p = ei{E) as derived in Section [3121 to a 
positive trigger flag. Taking all geometrical effects as 
well as the modular trigger efficiency into account, a 
detection efficiency of 

etot,MCveto-voIume = (93.6 it 1.5)% (5) 



is obtained. Out of the 6.4% of muons entering the 
veto volume undetected, more than a third of them 
(2.4% of all) pass through geometrical gaps of the 
system. This reflects the geometrical coverage of the 
//-veto system of 98%. In addition, there are 0.9% 
of all muons missed which pass through a module for 
which the data were not recorded during a small part 
of the total measurement time. This missing data is 
accounted for in the simulation by setting the trigger 
efficiency of this module to zero for the appropriate 
amount of simulated time. The statistical and sys- 
tematic uncertainty of the detection efficiency was 
evaluated assuming a set of different threshold val- 
ues for the modular trigger efficiency. Simulating the 
complete muon flux as described in Section 13. H it is 
noticeable that the number of triggered events is fur- 
ther increased by ^ 25% due to muons passing out- 
side the veto volume but having secondaries hitting 
at least one module. 

Similarly to the detection efficiency 
etot.MCvcto- volume, we also derived a detection 
efficiency for a much smaller volume to be pen- 
etrated by muons, i.e. a sphere with a radius of 
1 m around the center of the cryostat was chosen. 
This second case is intended as a more adapted 
efficiency value in consideration of the WIMP search. 
Whereas the first volume includes muons passing 
through some corners of the experimental setup with 
gaps due to mounting structures, the central sphere 
requires muon tracks that may produce secondary 
particles close to or within the cryostat. For these 
muons we obtain a detection efficiency of 

etot.MCcentralsphere = (97.7 it 1.5)% (6) 

Since we rely on the same calibration and mea- 
sured trigger efficiencies of the indiviual modules, 
the estimated error is of the same size as for 
etot.MCvcto- volume- The remaining inefficiency for this 
sample is almost entirely caused by the trigger inef- 
ficiencies of the individual modules and could thus 
be reduced by lowering the effective thresholds. As 
most of the muon-induced background in the bolome- 
ters is due to muons passing near the cryostat (see 
Section and Fig. [10|), etot,MCcentraisphere reflects 
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better the veto efficiency for potential muon-induced 
background than etot.MCveto-volume- 

3.4 Muon detection efficiency derived 
from bolometer data 

The following approach avoids the uncertainties of 
the energy calibration and MC simulation and in- 
stead purely relies on data acquired with EDEL- 
WEISS itself. It is fully independent of the meth- 
ods described in Section 13.31 and leads to the most 
reliable results for the muon detection efficiency . In 
this method, muon candidates are selected in the Ge 
bolometers and the efficiency is probed by a simple 
coincidence requirement. Above 6 — 7MeV of energy 
deposition in the bolometers, the fraction of events 
from the ambient radioactivity quickly drops and 
should be completely zero above 10 MeV for anything 
besides throughgoing muons. To ensure suppres- 
sion of all potential background in a single bolome- 
ter, an additional multiplicity requirement of at least 
two detectors was introduced. With a selection of 
-E'heat > 7 MeV and a multiplicity of mboio > 2, a to- 
tal of 34 ^ candidate events was found in the bolome- 
ters. All of these events were detected in the /i-veto 
system with time differences At = tboio — ^veto as can 
be seen in Fig. [71 As will be discussed in detail in 
Section 14. H positive values of At of the order of tens 
of ^s are expected for perfect synchronization of the 
bolometer and fi-veto system. Such differences are 
due to the sampling rate of lOOkS/s of the bolometer 
signals and the rise time of the ionization signal of 
the order of a few ^s. Larger differences can occur 
for high energy events which saturate the ionization 
channel and make the time reconstruction more dif- 
ficult. In contrast, there are 7 of the coincidences 
with an unphysical timing At < where the fast 
scintillation signal of the //-veto was observed up to 
12.1 ms after the signal in the bolometers (see Fig. [7] 
and Section 14.11 for the discussion of the synchroniza- 
tion accuracy). Most of the events including the 7 
events with unphysical timing, have reconstructed 
muon tracks (Section 13. 5p through the cryostat. A 
few of them show multiple hits indicating a shower 
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Figure 7: Time difference between bolometer muon candidates 
and their closest /i-veto event. Note that two events with At = 
—4.1 ms and —12.1 ms are not shown. 



in the whole system. Thus, with an expectation of 
0.14 events due to accidental coincidences in the time 
interval [—15, -|-1] ms, all 34 events are considered as 
muon-induced coincidences. 

The best estimate for the efficiency observing 34 
out of 34 events is etot,Data = 100%. We deter- 
mine a lower limit for the efficiency value by ex- 
tracting the lowest efficiency value for which one 
expects to see 34/34 events in 10% of the cases: 
P{e = x|34/34evts) = 10%. According to bino- 
mial statistics P{k,n,p) = (^)p'^(l — where p 
corresponds to the efficiency e and n = k = 34, we 
transform this into 



etot.Data ^ vO.l = 0.935 at 



C.L. 



(7) 



Thus, this method yields a lower efficiency limit 
of 93.5% entirely due to the statistical uncertainty. 
From the point of systematics, the candidate selec- 
tion is biased to close muon tracks. However, this is 
also the case for the dominant fraction of all muon- 
induced events in the bolometers as will be discussed 
in Sec l4.2[ In conclusion, the method described here 
is supposed to be adequate for determining the re- 
jection capability of the //-veto system. The sys- 
tematic uncertainty is estimated to be less than 1% 
and the lower boundary of 93.5% for the /i-veto sys- 
tem efficiency is regarded reliable albeit very 
weak boundary limited by the low statistics of the se- 
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lected bolometer event sample. It is in full agreement 
with the detection efficiency derived from MC as de- 
scribed in Section 13.31 With the further data taking 
of EDELWEISS, the statistics will increase and thus 
this method will also yield a more precise determina- 
tion of etot,Data- 

For all following investigations of muon-induced 
background rates, we will use the muon detection ef- 
ficiency derived here, i.e. etot = etot,Data- 

3.5 Muon track reconstruction 

The information provided by the module number and 
time difference of the two PMT signals on each end of 
a plastic scintillator module are used to reconstruct 
the coordinates of the energy deposition in the /i-veto 
system. Since there is no information about the po- 
sition of the event with respect to the module width 
of 65 cm and to the thickness of 5 cm, the position 
perpendicular to the module axis is set to a ran- 
dom value between and 65 cm in order not to cre- 
ate discrete steps in geometrical distributions of the 
recorded events. The thickness of the module is ne- 
glected. For each event generating coordinate-signals 
in two modules, a track can be defined by connecting 
the two penetration points, including their relative 
timing. This track allows to reconstruct azimuth and 
zenith angles as well as the minimum distance of the 
muon track to the center of the cryostat. 

The selection of events with exactly two energy 
depositions on different sides of the ^-veto system 
is supposed to be due to muons, since only those 
are penetrating enough to go through the whole ex- 
periment. The resulting angular dependent spectra 
are a negative imprint of the shielding properties of 
the mountain profile convoluted with the geometri- 
cal coverage and efficiencies of the experiment. In 
Fig. [8l the measured zenith distribution of muons 
traversing the EDELWEISS experiment is compared 
to MC simulations [2^1 [2S] and previous measure- 
ments of the Frejus experiment |41] . The MC simu- 
lation incorporates both the geometry and the single 
module efficiencies of the EDELWEISS experiment as 
determined in section 13.21 For better comparison of 
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Figure 8: Zenith angle distribution of muon tracks in the 
EDELWEISS fi-veto system (black data points), in the Frejus 
proton decay detector (dotted blue) [41] and of the correspond- 
ing simulated tracks (EDELWEISS solid red; Frejus dashed 
green). All spectra are normalized to the EDELWEISS data 
set. The simulation incorporates the full geometrical accep- 
tance functions. 

the angular distributions, all spectra are normalized 
to the EDELWEISS data set. The MC simulation 
reproduces well the EDELWEISS data (black data 
points). In contrast, the Frejus experiment [T7] has 
a different geometry and hence the acceptance is sig- 
nificantly different leading to a different zenith angle 
distribution as confirmed by simulated tracks with 
the according acceptance function. From Fig. [HI one 
can see that the differences in the measurements of 
EDELWEISS and Frejus are due only to the differ- 
ent acceptance functions. The simulated zenith angle 
distributions based on the same muon flux generator 
reproduce both angular distribution once the accep- 
tance of the experiments is incorporated. 

The actual profile of the EDELWEISS data with a 
maximum at 35° is well understood with a large rock 
overburden of 1800 m at zero degree and the small- 
est rock overburden at the fianks of the mountain. 
The spectral shape of the MC simulation is sensitive 
to the implemented energy thresholds of the individ- 
ual modules. The consistency of the simulation with 
the EDELWEISS data underlines that both the cal- 
ibration and single module efficiency determination 
as well as the track reconstruction with an estimated 
reconstruction uncertainty of fx^ = 15 cm along the 
module axis are well controlled. The zenith spec- 
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Figure 9: Azimuth angle distribution of muon tracks in the 
EDELWEISS /J- veto system (black data points), in the Frejus 
proton decay detector (dotted blue) [41] and of simulated tracks 
through EDELWEISS (solid red) and Frejus (dashed green). 
The Frejus detector has vanishing detection efficiency at and 
180 degrees. Frejus data and simulation results were normal- 
ized to the EDELWEISS data set. 
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Figure 10: Muon track distance to the center of the cryostat. 
EDELWEISS /i-candidates are shown as black data points and 
the MC simulation in solid red. 



trum has been cut at 80° since both the map of the 
mountain profile as well as the initial muon spectrum 
of the muon generator used in the MC have large un- 
certainties at such horizontal zenith angles. 

The azimuth angle distribution (Fig. [9|) also shows 
good agreement between the measured data in 
EDELWEISS and the MC simulation modehng de- 
tected muon tracks. The zero degree direction is 
aligned with the laboratory direction. In the geo- 
graphical coordinate system this corresponds to an 
offset of 16° to the north direction. There are two 
artificial peaks at 90° and 270° in the data. These 
peaks result from the fact that the modular position 
of events which have a TDC signal that would place 
them outside of a plastic scintillator module is set 
to be at the end of the module. For the top mod- 
ules, the module ends are in plane with the North 
and South vertical scintillator modules giving recon- 
structed tracks with 90° or 270° azimuth angle. The 
good agreement between the data and the simulation 
suggests that it describes well the spatial resolution 
and the proportion of tracks reconstructed near the 
module edges. 

The profile measured by Frejus is rather similar 
to this work with a few exceptions. One immedi- 
ately notes two dips in the angular spectrum. For 



azimuthal angles of 0° and 180° muons passed parallel 
to the detector planes of the Frejus experiment and 
could not be detected [17]. Again, the same muon 
generator reproduces the two experimental azimuthal 
spectra once the according acceptance functions are 
incorporated into the simulations. 

Another application of the /U track reconstruction 
is the identification of muons passing close to or pen- 
etrating the bolometers. The distribution of the min- 
imum track distance to the center of the cryostat 
is shown in Fig. [TOl The characteristic shape of 
the distribution originates from the number of muon 
tracks in a thin spherical shell with a certain fixed 
radius around the center of the cryostat and the fact 
whether this shell is contained within the /i-veto sys- 
tem or not. As the cryostat is 1.9 m away from the 
top and the four upper sidewalls, this distance corre- 
sponds to the maximum of the muon track distribu- 
tion. Afterwards the number of muon tracks within 
the spherical shell still increases, but the number of 
detected traces decreases again. The maximal dis- 
tance of a muon veto module to the center of the cryo- 
stat is (i ~ 4 m. The simulation reproduces very well 
the data and underlines the reliability of the muon 
track reconstruction. 
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3.6 Determination of the muon flux 

To conclude, we determine the integral muon flux 
through a horizontal surface. Requiring a module be- 
ing hit in at least two surfaces of the ^-veto system, 
we accumulated 18503 muon candidates which cor- 
responds to a measured rate of r^_cand. = (108.7 ± 
0.8)/d. Due to the rather cubic geometry and the 
individual detection efficiencies of each module, the 
overall detection acceptance for a given muon selec- 
tion must be determined through simulations. For 
the applied muon selection cuts, a combined detec- 
tion acceptance omc = 20.0 it 0.4candidates/(/i/m^) 
was derived from simulations accounting for the ge- 
ometry of the setup and modular efficiencies e^. This 
leads to an acceptance-corrected total flux of muon- 
induced events of 

^fi = ^ fi-cand. / auC 

= 5.4 ± 0.2(stat.)l°J(syst.) fi/m'^/d. (8) 

The systematic uncertainty is made up of three con- 
tributions. The first one arises from background 
events not induced by muons but being in coincidence 
in modules of two geometrical surfaces. This may 
happen in corners of the //-veto system. To check 
their influence, we analyzed subsamples where adja- 
cent surfaces are not allowed. This reduces r^_(,and. 
as well as omc leading to an overall reduction of 
by 10%. As described in Section [3.21 the modu- 
lar efficiencies depend on the calibration param- 
eters. Allowing a variation of ±20 % leads to an- 
other 10 % uncertainty. Third, there is a remain- 
ing uncertainty of the derived flux due to modeling 
the muon flux based on single muons produced on 
the hemisphere and neglecting muon bundles which 
were measured in the Frejus experiment contribut- 
ing 5% to the overall muon flux [17] . In our sim- 
ulations, the average number of muons per event is 
only slightly larger than one due to secondary muons 
produced in the rock between the start hemisphere 
and the lab. Deriving alternatively from the 18503 
measured muon candidates the muon flux through a 
sphere leads to = 16.4 it 0.4 candidates/ (/i/m^) 
and subsequently ^^iF^ = Q.Q ^/m? / d. 



The measured flux as given in equ. [8] is slightly 
higher but yet compatible with the measurements 
of the Frejus experiment. To compare the values, 
the different definitions and cuts have to be taken 
into account. The event rate determined here cor- 
responds to the sum of events with muon multiplic- 
ity m > 1, including muon bundles. In addition, 
no zenith angle cut of < 60° as in [T7] was ap- 
plied in our measurement. Since it is not stated 
in |17) whether the muon flux was measured for a 
horizontal area or a sphere, we give the appropri- 
ate values for both assumptions here. Interpreting 
the quoted flux in [17| flux through a horizon- 
tal plane and correcting it for the inefficiency of the 
60° cut, evaluated from our measurements, a value of 

^F?ljur*^^ = (5.2±0.1stat.)Ai/mV(i is calculated. The 
alternative interpretation of the values in |17] as the 
flux through a sphere would lead to a smaller flux of 

^F?6jur*''^ = (4.2 ± 0.1stat.)M/mVc^ to be compared 
with the value given in equ. El 



4 Coincidences in the /i-veto sys- 
tem and bolometers 

For WIMP searches, every background component 
has to be well understood and maximally suppressed. 
In this section, the muon-induced event population in 
the bolometers is studied in situ, for the specific geo- 
metrical arrangement of bolometers during the data 
acquisition period from March 2009 to May 2010. 
In addition to quantifying the muon-induced neutron 
component for the WIMP analysis [1^, the topology 
of muon-induced events is investigated. With this 
derived information, an estimate of muon-induced 
neutrons can be deduced for the next phase of the 
experiment, EDELWEISS-III. The investigation is 
performed using a newly developed data framework, 
KData |22) that allows the event-based storage and 
analysis of both bolometer and /i-veto data in the 
same structure. 
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Figure 11: Time diflerence between bolometer events and hits in the /i-veto system, (a) data taken after September 2009, (b) 
data taken from March to September 2009. Black indicates coincidences while the blue dotted line corresponds to the expected 
rate of accidental background. The vertical lines represent timing cuts applied for further investigation of coincidences. Solid 
lines indicate the applied cut for the appropriate data period while the dotted lines describe the cut used for the other data 
period. 



4.1 Determination of the coincidence in- 
terval 

Searching for delayed coincidences between a prompt 
muon veto hit and a secondary bolometer event, a 
typical time delay At = tboio — ^veto of a few tens 
of microseconds is expected. This time difference is 
mainly due to the rise time of the bolometer pulses 
and their sampling rate of lOOkS/s. The accuracy of 
the synchronization of the separate acquisition sys- 
tems leads to a broadening of a few tens of microsec- 
onds. The overall recorded time difference as shown 
in Fig. [11] is thus not caused by physics processes 
such as thermalization of neutrons. The exact exper- 
imental coincidence window to select muon-induced 
signals is determined by the At distribution itself. 
In the following, we select events where at least two 
modules were hit (rriveto ^ 2 with at least one ADC 
or TDC entry per module being non-zero) leading to 
an event rate of 0.12 Hz. This is a much stricter re- 
quirement than in vetoing potential backgrounds in 
the dark matter search where bolometer events are 
rejected whenever there is any event recorded in the 



/i-veto system within a it 1ms time interval. With 
an event rate of 0.41 Hz in the /i-veto system, the 
applied dead time amounts to less than 0.1 % of the 
total measuring time. Having a muon candidate with 
"T-voto > 2, we then extract the time difference At to 
the nearest event in the bolometer data set. 

This time difference is shown in Fig. [TT] for two 
distinct periods of data taking. Fig. [TTk shows the 
distribution for the data after September 2009 with 
fully functioning synchronization between /i-veto sys- 
tem and bolometers. A clear peak of muon-induced 
coincidences can be identified and a time interval 
of 20 /is < At < 80 fis was chosen maximizing the 
signal-to-background ratio containing at least 97.5 % 
of the peak population. For the early data set from 
March 2009 to mid-September 2009 corresponding to 
about 1 /3 of the total exposure, some data losses in 
the optical transfer of the 48-bit clock from bolome- 
ters to the /i-veto system degraded the synchroniza- 
tion performance. An offline correction to the syn- 
chronization could reproduce a peak in At, but at 
slightly unphysical time differences (Fig. Illb). The 
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corresponding interval for coincidences was derived to 
be -300 fis < At < O^s. However, not all of the syn- 
chronization errors could be corrected in the offline 
analysis. A detailed investigation of coincidences 
with synchronization problems resulted in the quan- 
tification of (12 lb 8) % of all cases with At < —1ms 
(including the two outliers at At = —12.1ms and 
At = —4.1ms in Section [3^ . An attributed un- 
certainty of a few ms in synchronizing the systems 
is much larger than the vetoing interval of it 1ms 
around each veto hit. When deriving the muon- 
induced background in the WIMP search in Sec- 
tion 231 we therefore conservatively treat these peri- 
ods with an exposure of (17 ± 11) kg • d as bolometer 
data without functioning /i-veto system. 

With the two coincidence intervals as shown in 
Fig. nil a total of 76+158 = 234 events summing over 
the 1st and 2nd period were selected for further in- 
vestigation. An accidental background of 6.75 it 0.03 
events was determined in a two second wide window 
with an offset of 5 seconds before the veto event. 
This corresponds to an achieved signal to noise ratio 
of better than 33/1. After subtraction of accidental 
background, the remaining excess of coincidences in 
EDELWEISS-II is entirely due to events induced by 
muons. 

4.2 Topology of muon-induced bolometer 
events 

A WIMP scattering off a Ge nucleus would lead to a 
single scatter event in a bolometer, i.e. a bolometer 
multiplicity rn-boio = 1- Any event with rn-boio > 1 can 
therefore be rejected. The multiplicity of bolometer 
hits in a bolometric array is thus a potentially pow- 
erful tool to reject background, depending on the na- 
ture of the background as well as on the total detec- 
tor mass and the granularity of the crystal array. The 
multiplicity of the 234 extracted /i-induced bolometer 
events is shown in Fig. [T2[ For the specific setup of 
EDELWEISS-II, (46.6 ± 0.6) % of all muon-induced 
events in the complete energy range would be rejected 
in a WIMP search due to mboio,^ > Ij with an aver- 
age bolometer multiplicity of (mboio.^i) = 2.27 + 0.18. 
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Figure 12: Bolometer multiplicity for muon-induced events 
(solid red). For comparison, a scaled histogram of the multi- 
plicity of all bolometer events is superimposed (dotted black). 

This is in clear contrast to (?Tiboio,aii) = 1-05 for all 
bolometer events, with the rriboio.aii distribution also 
shown in Fig. [T2j 

The fraction of events with mboio.^t > 1 strongly 
depends on the geometrical setup. This can be seen 
in a comparison of the present results with the sim- 
ulations based on a fully equipped cryostat with 120 
Ge crystals, each with a mass of 320 g [26J, adding up 
to 38.4 kg. As described in Section l2.ll the current 
study incorporates only 10 ID detectors with a total 
mass of 3.6 kg. In the ideally packed array used in the 
simulations, a fraction of 84% of all muon-induced 
events have a multiplicity of mboio,^ > 1- Thus, with 
increasing numbers of bolometers in a densely packed 
array, a muon-induced background can be efficiently 
suppressed due to its high multiplicity. 

Furthermore, muon-induced background events ex- 
hibit an energy spectrum which is very different to 
an exponentially decreasing recoil spectrum expected 
from WIMP scattering. Fig. [13] shows the sum en- 
ergy of all energy depositions in the detectors asso- 
ciated with a muon-induced event. The spectrum is 
composed of two main regions. Up to 1 MeV the en- 
ergy deposition in the bolometers is caused by single 
or multiple energy depositions of secondaries, mainly 
electrons and gammas. From 50keV to several 100 
keV the dominant interaction is Compton scattering. 
The second feature is a broad peak between 5 and 
100 MeV with a maximum around 20 MeV caused by 
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Figure 13: Sum of energy deposition in the bolometers for 
the 234 coincidence events (dashed black - 4 bins per decade) 
as well as for muon-induced single bolometer hits (long-dashed 
green) and background events in individual bolometers (dotted 
blue - 64 bins per decade). For comparison, a subsample of 109 
events is shown for which a muon track could be reconstructed 
(solid red, see text for more details). 

muons passing through the Ge crystals (energy loss 
of muons in Ge ~10MeV/cm). For those events, 
the energy deposition is mainly dependent on the 
path length of the muon in the crystals. For sin- 
gle scatter events (Fig. [13] long-dashed green) this 
feature is strongly suppressed, indicating that these 
are mainly caused by secondaries. Last, there is 
a small fraction of events below 50keV, where the 
dominant component is elastic scattering of muon- 
induced neutrons as will be discussed later. Fig. [13] 
also shows the energy deposition in the bolometers 
for background events (dotted blue histogram) which 
clearly contrasts the muon-induced distribution. The 
peaks in this spectrum can be associated to a decays 
of ^^°Po (5.4 MeV), which can be further suppressed 
by a fiducial volume cut, a Compton backscattering 
bump above 100 keV as well as the lOkeV gamma 
line from cosmogenic activation. 

Apart from the bolometer multiplicity and energy 
deposition of muon-induced events, the distance of 
the muon track from the center of the cryostat could 
be deduced from the information in the muon veto 
modules. Fig. [TJ] shows the distance of 109 out of 
the 234 sequential events. For these 109 events, the 
muon track could be reconstructed unambiguously 
due to full spatial information in exactly two of the 
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Figure 14: Distance to the center of the cryostat for 109 co- 
incidence events with track reconstruction (solid black), with 
energy depositions in the bolometers below (dotted blue) and 
above 1 MeV (dashed red). Also shown is the scaled distribu- 
tion of muon events identified in the /i-veto system without 
requiring a coincidence (dash-dotted green) with its maximum 
at f» 2 m as described in Section 13.51 and Fig. 1101 

veto modules. More than 90% of these muon-induced 
bolometer events are produced by muons passing 
within a distance of less than Im to the cryostat's 
center whereas the overall distribution of all identi- 
fied muons is much broader with a maximum around 
d = 2 m. Although muon-induced bolometer events 
have short track distances to the cryostat, it is impor- 
tant to note that this does not translate directly into 
muons hitting the Ge crystals with large energy de- 
positions. As can be seen in Fig. [TJ] the distances for 
events with energy depositions in the bolometers be- 
low (dotted blue) and above 1 MeV ( dashed red) 
have almost identical distributions. Thus, muons 
even passing through the cryostat lead to energy de- 
positions below 1 MeV, reaching the typical WIMP 
range of tens of keV, indicating energy deposition by 
secondary particles. On the other hand, distances of 
d > 0.7 m correspond to muons passing outside the 
cryostat, where secondary electrons and gammas are 
unlikely to reach the Ge bolometers due to the var- 
ious shields. In this case, coincident events can be 
interpreted as muon-induced neutrons scattering off 
the Ge nuclei. 

The requirement for reconstructing muon tracks 
may introduce a slight bias in selecting muon-induced 
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Figure 15: Simulation of production vertices of neutrons in 
muon-induced showers leading to an energy deposition by nu- 
clear recoils in the bolometers [2^. Shown is a side projection 
of the underground laboratory LSM. Most of the neutrons are 
produced in the Pb shield and the steel of the support struc- 
ture. 

events with high energy deposition in the bolometers 
as can be seen in Fig. [T3l Low energy depositions 
originating from muon-induced neutrons might also 
be connected to muons passing outside the /i-veto 
system. However, the production vertices for these 
neutrons are dominantly positioned in the Pb shield- 
ing as can be seen in Fig. [T5j From these simulations, 
one expects only 0.05 % of all neutrons produced out- 
side of the /i-veto system, while more than 90 % are 
produced within the lead shield or the cryostat it- 
self. In conclusion, muon-induced bolometer events 
(including potential neutron scattering events) are 
mainly caused by muons with a track of less than 
1 m distance from the center of the cryostat and thus 
potentially can be vetoed with the //-veto system. 

4.3 Muon-induced background for dark 
matter search 

The muon-induced background for the WIMP analy- 
sis is determined as the product of the rate of muon- 



induced single hit events in the nuclear recoil band, 
the accumulated exposure of 384 kg • d and the fi- 
veto system inefficiency (1 — etot) as determined in 
eq. d?]). For a DM search, it is crucial to have the 
highest muon detection efficiency possible. Hence, in 
this analysis we abandon the requirement for muon 
candidates as used in Sections |4. 1 1 and i.e. having 
hits in more than one veto module. The total num- 
ber of coincidences then moderately increased from 
Ar™>i = 234 to Nq = 283 events, while the rate of 
accidental background increased by a factor of 3.5 
from N^>^ = 6.75 ± 0.03 to Na = 23.61 ± 0.06. As 
for the 234 coincidences studied before, there was no 
rejection of periods with high baseline noise nor any 
selection of fiducial bolometer events. Thus, the total 
exposure is much larger than in the WIMP analysis 
and amounts to 1504 kg • d. With a net number of 
N = Nc — Na = 259 it 17 coincidences, the rate of 
muon-induced events in the complete energy range, 
not discriminating between electron and nuclear re- 
coils, amounts to 

pM-ind. cvts ^ (Q^^72 ± 0.001) evts/(kg • d). (9) 

The rate of muon-induced background events within 
EDELWEISS-II is in good agreement with earlier in- 
vestigations [21] based on a much smaller event sam- 
ple and slightly different selection cuts. Fig. [16] shows 
the events as individual bolometer hits with an en- 
ergy deposition below 250 keV. Applying a fiducial 
volume cut inferring reliable charge collection and 
thus ionization yield and a x^"Cut comparing heat 
and ionization pulse fits with their template pulses, 
bolometer hits shown as dots are rejected. After a 
time-averaged 90 % C.L. nuclear recoil band selection 
applied for each single detector, seven bolometer hits 
remain (see Fig. [16] within the red bands). Out of 
these seven, six are above the 20 keV analysis thresh- 
old, and only four (shown as red crosses) are single 
bolometer hits without a coincident hit in another 
bolometer as required for a WIMP candidate. (The 
bolometer hit at 60 keV is in coincidence with another 
one, the hit at 70keV is outside the predefined nu- 
clear recoil band.) Three of these events have recon- 
structed muon tracks going through the lead shield 
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Table 1: Rates of coincidences. All results are given after back- 
ground subtraction. The cuts on recoil energy E^e* and Q- value 
Qrec are fulfilled, when at least one bolometer hit within an 
event fulfills the condition. For fiducial nuclear recoil (NR) 
events, the 90 % signal region is defined following the Lindhard 
parametrization 1331. 



Selection CLit r^_ind (cvts/kg/d) Exposure (kgd) 



None 0.172 ±0.012 1504 

E^^'< 250keV 0.095 ± 0.009 1504 

E^^;; < 250 keV; Q"* < 0.6 0.035 ± 0.005 1504 

eI*^' < 250keV; fid. NR 0.017 ± 0.006 481 

E''^' < 250 keV; fid. NR; single O.OOStpoM 481 



Figure 16: Scatter plot of bolometer hits in the coincidence 
analysis for energies below 250 keV. Fiducial events that are ex- 
pected to have a full charge collection and thus a reliable value 
for Q are drawn as crosses (blue and red) . The events marked in 
red satisfy all WIMP criteria. Other events are shown as black 
dots. The lines correspond to the 90 % C.L. gamma region 
(blue), the 90% C.L. nuclear recoil band (red), the analysis 
threshold of 20 keV used in the WIMP analysis (black) and an 
ionization threshold for fiducial events of 3 keV (green) . Except 
for the analysis threshold, all resolutions and bands are time 
and bolometer dependent. For illustration, the time-averaged 
bands of an exemplary detector, ID401, are drawn here. 

and even through the cryostat itself. The fourth 
event is in coincidence with an event in the /i-veto 
system, where 11 plastic scintillator modules have 
been hit, indicating a muon-induced shower. Fig. [16] 
shows individual bolometer hits. Rates in Table [T] 
refer to bolometer events with potentially more than 
one bolometer hit. 

Applying all cuts required for WIMP candidates 
results in an effective fiducial exposure of 481 kg-d. 
Thus, the rate of muon-induced WIMP-like events 
amounts to 

r^- = (0.0081°:°°^) evts/(kg-d). (10) 



calculated according to 

- A/^+^r'^-"(l - etot) + Mlxpr'^"" = 0.40 (11) 

Here M^^p = 38 ± 1 1 kg • d denotes the sum of exposure 
of a short period where only bolometer data was acquired 
(21 kg-d) and the earlier period where a hardware mal- 
function led to a faulty synchronization in 12 % of the 
measuring time (17±llkg-d). M^+^ = 384 kg • d - M|^p 
is the exposure of the WIMP search data set, when both 
acquisition systems (bolometer and /z-veto system) were 
running with full synchronization. The above expectation 
value is dominated by the contribution from M^^^. 

Taking into account all uncertainties on F**^", etot and 
Afg^p we derive an upper limit of iV''^" < 0.72 events 
(90%C.L.). Comparing this result with other back- 
ground estimates for the EDELWEISS-II WIMP search 
from misidentified gamma events N'^ < 0.9, from surface 
events < 0.3 [13] and from ambient neutrons < 3.1 
[14] , muon-induced background contributes less than 20 % 
to the total background. Note that this is dominantly due 
to including periods in the WIMP search where the /i-veto 
system was not synchronized properly to the bolometer 
system. For future searches as in EDELWEISS-III, such 
periods will not be included and thus an even better re- 
jection of muon-induced background can be achieved. 



Assuming a conservative veto efficiency as in Eq. \7\ 
the rate of muon-induced WIMP-like events amounts 
to an irreducible background rate of P^^""^ < 5.8 • 
10-4 evts/(kg • d) at 90 % C.L. 

From Equations [7] and [TO] for the fi-veto detection effi- 
ciency and the rate of muon-induced WIMP-like events, 
the number of expected unvetoed muon-induced single 
scatter neutron events A^^~" within the WIMP analysis is 



5 Expected background for 
EDELWEISS-III and EURECA 

For the forthcoming phase of EDELWEISS-III as well as 
future experiments, an accurate number for the muon- 
induced neutron background can be predicted only with 
full modeling of the specific geometry due to the suppres- 
sion of multiple scattering signatures in a densely packed 
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detector array and the benefits from additional sliielding. 
A larger number of detectors, such as the forty 800 g detec- 
tors planned for EDELWEISS-III, will increase the gran- 
ularity and hence the number of multiple bolometer hits. 
Also, the effective area of Ge bolometers exposed to the 
muon flux remains approximately constant while the de- 
tector mass increases proportionally with the number of 
detector planes. Furthermore, the introduction of an ad- 
ditional PE shield inside the cryostat is going to moderate 
the neutron flux induced in the lead. All effects will reduce 
the rate of muon-induced WIMP-like background events. 

Taking the measured rate of muon-induced nuclear re- 
coil events (Eq. [TU|) and the derived efficiency of the /i-veto 
system of 97.7 % as in Eq. M a total of iVg^ooofcg d = ^-^-al 
(90%C.L.) irreducible background events due to muon- 
induced neutrons are projected for EDELWEISS-III for 
a 6-month exposure of 3000 kg • d. Both the background 
rate and the muon detection efficiency will be measured in 
the ongoing analysis of the EDELWEISS-III experiment 
and in parallel modeled with its specific new setup and 
better shielding. 

Going further to ton-scale experiments such as the en- 
visaged EURECA experiment, the muon-induced neutron 
background could be of sizable contribution. The strategy 
of reducing it with respect to the EDELWEISS design is 
to reduce the high-Z material close to the bolometers as 
well as including a highly efficient water Cerenkov active 
shield around the cryostat [5T| . 

6 Conclusion and outlook 

The angular dependence of the muon flux through the 
EDELWEISS experiment has been analyzed and is well 
modeled by MC simulations. Two independent methods 
to determine the /i-veto system efficiency have been pre- 
sented, giving results limited so far by small statistics. 
A conservative value of etot > 93.5% has been used for 
the calculation of muon-induced background in the WIMP 
search data set. The data of the bolometer array together 
with the modular structure of the /i-veto system allowed 
a detailed study of the topology of muon-induced events 
within EDELWEISS. The analysis of muon-induced back- 
ground has been performed in the specific setup of the data 
acquisition of 2009/2010 with 10 Ge bolometers. For this 
configuration, the rate of muon-induced events that mimic 
WIMPs has been determined to F'^-" = (0.008i|^:EJ|^^) 
events/ (kg • d). In conclusion, the muon-induced back- 
ground with an expectation of N'^~" < 0.72 events within 
this WIMP search was small compared to other back- 



grounds. Thus, only a minor fraction of the five events 
observed in the dark matter search [13 is assumed to be 
due to muon-induced interactions. 

The feasibility of a background free data taking of at 
least 3000 kg -d (EDELWEISS-III 6 months) has been 
shown with a conservative extrapolation from current 
data. Expected benefits from larger bolometer granularity 
and hence larger multiplicity of events as well as from the 
added shielding have been briefly discussed and should 
enable background free operation beyond a sixth month 
exposure. However, future ton-scale experiments, like the 
EURECA experiment, need different shielding concepts 
with less high-Z material close to the bolometers. 
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